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1 ^ ■ phase and the superfluid gap in the low-density regime are computed, using a realistic Hamiltonian 
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[2 ■ INTRODUCTION 

O 

VO \ The improved accuracy of experimental data on nuclei, together with a rediscovered role 

t~^ ■ of nuclear matter properties in the understanding of nuclear structure and several phe- 

O \ nomena of astrophy sical interest[ 1 ] , has shown the need for a more detailed investigation 

(^ \ of the ground state of nuclear many-body systems. 

Recent realistic nuclear Hamiltonians have been used to compute properties of light 

. nuclei in very good agreement with experiments [2]. However, the physical properties 

rN \ of nuclear and neutron matter could be very different from that of nuclei; in fact the 

c^ ■ nucleon density in the core of heavy nuclei reaches the maximum value of po=0.16 

fm^^, while the relevant range of density of matter inside neutron stars is up to 9 times 

PolbO. Therefore we are now facing, on one side, the problem of finding a fundamental 

scheme for the description of nuclear forces, valid from the deuteron up to dense nuclear 

matter, which is still an open fundamental problem, and, on the other, that of solving 

a many-body problem which is made extremely complex by the strong spin-isospin 

dependence of the forces. In this paper we will not address the problem of determining 

the nuclear force. We will consider a nuclear Hamiltonian which provides good fits to 

the N-N data up to meson production and reproduces fairly well the ground state and 

the low energy spectra of light nuclei. This is made of two- and three-body spin-isospin 

dependent potential. We present results for the ground state of large nuclear systems 

with this Hamiltonian. 

It is well known that Quantum Monte Carlo (QMC) methods can provide estimates 



of physical observables at the best known accuracy[4], and they are therefore useful to 
gauge the validity of proposed interaction models without having the bias of using more 
approximate methods. 

A new generation of powerful QMC techniques have been recently devised to simu- 
late large nucleonic systems with up to hundred of nucleons: the Auxiliary Field Diffu- 
sion Monte Carlo (AFDMC)[5]. They have been used to compute the EOS of nuclear 
matter[|6|], showing important limitations of other many-body methods and of the mod- 
em nuclear interactions based on two- plus three-body potentials in the high density 
regime. The accuracy of AFDMC was demonstrated by comparing the ground state of 
light nuclei with results provided by Green's Function Monte Carlo (GFMQdVi] that is 
known to give accurate results for the properties of light nuclei up to A=12|l2|]. By sam- 
pling the spin-isospin states of the nucleons AFDMC can be applied to large systems; 
it was used to simulate the ground state of medium sized nuclei up to A=40[8], nu- 
clear matter with up to A=108[6], the properties of neutron-rich nuclei[i9|,il^, neutron 
drops ini], and neutron matter with up to A=l 14[[Tll[l2l.[T3|]. 



In this paper we discuss the latest results of the computation of the equation of state 
(EOS) of neutron matter in both the normal and superfluid phases, the computation of 
the superfluid gap of neutron matter in the low-density regime, and some preliminary 
results about the EOS of isospin-asymmetric nuclear matter. 

HAMILTONIAN 

The ground state of nuclear systems can be realistically studied by starting from the 
non-relativistic nuclear Hamiltonian of the form 

where m is the averaged mass of proton and neutron, and vtj and Vijk are two- and 
three-body potentials; it seems that the effect of forces due to n-body terms with 
n > 3 in the low energy properties of light nuclei is negligible. Such a form for the 
Hamiltonian has been shown to describe properties of light nuclei in excellent agreement 
with experimental data (see Ref. [2] and references therein). All the degrees of freedom 
responsible for the interaction between nucleons (such the ;r, p. A, etc.) are integrated 
out and included in v,j and Vijk- 

At present, several realistic nucleon-nucleon interactions (NN) fit scattering data with 
very high precision. We consider the NN potentials belonging to the Argonne family. 
Such interactions are written as a sum of operators: 



M 



\0^'\iJ), (2) 



where O^P'{i,i) are spin-isospin dependent operators. The number of operators M 
characterizes the interaction; the most accurate of them is the Argonne AVI 8 with 
M=18 [ | 14j] . Here we consider some simpler forms derived from AVI 8, namely the AV8' 



and the AV6' [15] with a smaller number of operators. For many systems, the difference 
between these simpler forms and the full AVI 8 potential can be computed perturbatively. 
Most of the contribution of the NN is due to the one pion exchange between nucleons, 
but the effect of other mesons exchange as well as some phenomenological terms are 
included. 

The eight O^P\i,j) terms in AV8' are given by the four central components 1, f,- ■ fj, 
di ■ Oj, {Oi ■ Oj){Xi ■ Tj), the tensor Stj, the tensor-T component 5,/f/ ■ fj, where 5// = 
3{di-rij){dj-rij) —di-dj, the spin-orbit L,j ■ 5,^ and the spin-orbit-T (Lij-Sij){Xi- Xj), 
where Lij and Sij are the total angular momentum and the total spin of the pair ij. 

The AV6' has the same structure of AV8', but the spin-orbit operators are dropped. 
In general, all the AVx' interactions are obtained starting from the AV18, written by 
dropping less important operators, and refitted in order to keep the most important 
features of NN in the scattering data[15]. 

The three-body interactions (TNI) is essential to overcome the underbinding of nuclei 
with more than two nucleons. The NN is fitted to scattering data and correctly gives the 
deuteron binding energy, but starting with ^H the NN is not sufficient to describe the 
ground state of light nuclei. The Urbana-IX (UIX) potential corrects this limitation of 
NN, and was fitted to light nuclei and to correctly reproduce the expected saturation 
energy of nuclear matter[16]. It essentially contains the Fujita-Miyazawa term[17] that 
describes the exchange of two pions between three nucleons, with the creation of an 
intermediate excited A state. Again, a phenomenological part is required to sum all the 
other neglected terms. The generic form of UIX is the following: 

V,jk = Vin + VR. (3) 



The Fujita-Miyazawa term 11 1711 is spin-isospin dependent: 



Vm = A2n £ {Xij,Xj^}{Xi ■ Xj, Xj ■ tJ + T [Xij,Xjk] [Xi ■ Xj, Xj ■ Xk] 



eye 



(4) 



where the Xjj operators describe the one pion exchange, and their structure is the same 
of that of AV6'. The phenomenological part is 

V,% = UoY^T\m^nj)T^{m^rjk) . (5) 

eye 

The factors A2k and Uq are kept as fitting parameters. The binding energy of symmet- 
rical nuclear matter is not well reproduced by such force. Other forms of TNI, called 



Illinois forces Hi 811 . which includes three-nucleon Feynman diagrams with two Deltas 
intermediate states, are available. However, they provide unrealistic overbinding of neu- 
tron systems when density increases lili [120 and they do not seem to describe realisti- 



cally high density (already at p > po) nucleonic systems. 

THE AFDMC METHOD 

Ground state AFDMC simulations rely, as do other traditional QMC methods, on pre- 
vious variational calculations, often performed within FHNC theory [ | 19|1 . to compute a 



trial wave function ^j-, which is used to guide the sampling of the random walk. A typ- 
ical form for ^7 is given by a correlation operator F operating on a mean field wave 
function 4>(i?), 

{R,S\^'t)=F^{R). (6) 

Mean field wave functions 4>(7?) that have been used are: (i) a Slater determinant ^fg 
of plane wave orbitals for nuclear and neutron matter in the normal phase, (ii) a linear 
combination ^^p of a small number of antisymmetric products of single particle orbitals 
^j{fi,Si) for nuclei and neutron drops, and (iii) a pfaffian ^p/, namely an antisymmetric 
product of independent pairs for neutron matter in superfluid phase. 

A realistic correlation operator is the one provided by FHNC/SOC theory, namely 
'^Wi>iYJ^=i f {nj)O^P' (?', j), where y is the symmetrizer and the operators O^^^ (/, j) 
are the same as those appearing in the two-body potential. 

Unfortunately, the evaluation of this wave function requires exponentially increasing 
computational time with the number of particles. This procedure is followed in varia- 
tional and Green's function Monte Carlo calculations, where the full sum over spin and 
isospin degrees of freedom is carried out. Since for large numbers of particles it is not 
computationally feasible to evaluate these trial functions, the much simpler correlation 
operator Y[j>if^{fij)^ which contains the central Jastrow correlation only, is used in- 
stead. The evaluation of the corresponding trial function requires order A^ operations to 
evaluate the Slater determinants and A^ operations for the central Jastrow. Since many 
important correlations are neglected in these simplified functions, we use the Hamilto- 
nian itself to define the spin sampling. 

The AFDMC method works much like Diffusion Monte CarloflSHIzilS]. The 



wave function is defined by a set of what we call walkers. Each walker is a set of the 3A 
coordinates of the particles plus a number A of four component spinors each representing 
a spin-isospin state. The imaginary time propagator for the kinetic energy and the spin- 
independent part of the potential is identical to that used in standard diffusion Monte 
Carlo. The new positions are sampled from a drifted Gaussian with a weight factor for 
branching given by the local energy of these components. Since they do not change the 
spin state, the spinors will be unchanged by these parts of the propagator. 

To sample the spinors we first use a Hubbard Stratonovich transformation to write the 
propagator as an integral over auxiliary fields of a separated product of single particle 
spin-isospin operators. We then sample the auxiliary field value, and the resulting 
sample independently changes each spinor for each particle in the sample, giving a new 
sampled walker. 



More details about the AFDMC method can be found in Ref. [[llfl 



NUCLEONIC MATTER 

The properties of nuclear matter, like the Equation of State (EOS), are of fundamental 
importance in nuclear physics, mainly because nuclei behave very much like liquid 
drops. Indeed, each of these can be associated with a mass formula, which fits the 
corresponding data of stable nuclei from A ~ 20 on. Any such mass formula has a 
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FIGURE 1. The AFDMC equation of state evaluated by simulating 66 neutrons in a periodic box. See 
the text for details. 



volume and a symmetry term provided by symmetrical nuclear matter and nuclear matter 
with N > Z respectively. Moreover, accurate model independent calculations of the 
above observables are much needed in the physics of heavy ion reactions, as well as 
in that of lepton and neutrino scattering off nuclei at intermediate energies. Medium 
effects have to be taken into account for the data analysis of such reactions at the present 
level of accuracy. 

In addition, the theoretical knowledge of the properties of asymmetric nuclear matter 
at low temperature is needed to predict the structure, the dynamics and the evolution of 
stars, in particular during their last stages, when they become ultra-dense neutron stars. 

We present in the following the results obtained with AFDMC for the EOS of pure 
neutron matter in normal phase, as well as the gap of the BCS phase of neutron 



matterflill |22L llj 
found in Refs. m 



Previously results of the EOS of nuclear matter and nuclei can be 

Neutron matter is simulated by considering A^^ neutrons in a periodic box, and particu- 
lar care is taken to evaluate the effects due to the finite size of the box. More details can 
be found in Refs. iH [H 

In Fig. [H we plot the AFDMC equation of state, obtained with the energy of 66 
neutrons, and the variational calculation of Akmal et al. of Ref. [23], where the AV18 
NN interaction combined with the Urbana UIX TNI was considered. As it can be seen 
both the AV8' and the AVI 8 essentially give an EOS with the same behavior, but the 
addition of the TNI adds some differences, in particular at higher densities. 

The AV8' interaction should be more attractive than AVI 8 as shown in light nuclei 
and in neutron drop calculation[7]. This result is not confirmed by our calculations as 
it is clearly visible in Fig. [T] The AFDMC has proved to be in very good agreement 
with the GFMC results for light nuclei[8], and we believe that the same accuracy 
is reached in the neutron matter calculation, as shown in the comparison with the 
GFMC results of 14 neutronsUli Usl]. On the other hand the AFDMC calculation of 
the nuclear matter have shown that FHNC/SOC does not seem to provide safe energy 
upperbounds. This is because of the lack of cluster diagrams with commutator terms 
beyond the SOC approximation and that of elementary diagrams, not included in the 
FHNC summation!^. 




FIGURE 2. The EOS of neutron matter in the low-density regime. The two calculations were performed 
using different trial wave functions modeling a normal and a BCS state. 



The addition of the UIX three-body interaction to the Hamiltonian increases the dif- 
ferences between the AFDMC results and that of Akmal et al. The difference cannot be 
due to finite-size effects in our calculation for the following reason: the total contribu- 
tion of UIX should be positive in neutron matter, so that the inclusion of box corrections 
as done for the two-body part of the Hamiltonian would eventually increase the to- 
tal energy. The periodic box-FHNC estimation of these effect essentially confirms this 
observation[12]. 

It is worth observing how important the three-nucleon interaction already is at 
medium-high densities. Its contribution at 2po is ~ 25MeV and increases very rapidly 
with density. The four Illinois potentials[18], built to include two A intermediate states 
in the three nucleon processes, lead to very different results compared to the Urbana 
IX III2I [Till EOS at medium-high densities, in spite of the fact that all of them provide 
a satisfactory fit to the ground state and the low energy spectrum of nuclei with A < 8. 
This, once more, points outs the importance of understanding the role of n-body forces 
with n > 3 in nuclear astrophysics. 

We explored the superfluid phase of low-density neutron matter. Because the AFDMC 
projects out the lowest energy state with the same symmetry and phase as the trial wave 
function, we tried to repeat some calculation using a BCS trial wave function of the 
form of Ref. ['24']. The BCS state provides an energy that is lower than the normal 
state, however the difference never exceeds 3% of the total energy. We report the energy 
per neutron in the low-density regime evaluated using a normal and a BCS trial wave 
function in Fig. [2l 

In the low-density regime the neutron matter is in a ^ So superfluid phase. In this 
regime, the neutron-neutron interaction is dominated by this channel, whose scattering 
length is very large and negative, about a=-18.5 fm. Several attempts to compute the 
pairing gap using a bare effective interaction like those used for cold atom problems have 
been performed in the last few years[22]. However, our results show that an accurate 
calculation of the superfluid gap in this regime must include the full Hamiltonian instead 
of one describing the ^Sq physics only. 

We computed the superfluid gap by means of AFDMC that allows for quantum 
simulations of the superfluid phase of neutron matter, by solving the ground state of 
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FIGURE 3. The superfluid gap computed using AFDMC and compared with other techniques is 
displayed. 
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FIGURE 4. The EOS of nuclear matter as a function of the density and of the isospin-asymmetry. See 
the text for details. 



the full Hamiltonian with AV8'+UIX without the use of some simplified interaction. In 
Fig. [3] the AFDMC gap is compared with standard BCS theory and some other often 
referred calculations based upon to correlated theories at two-body approximation (see 
also Ref. iH). 

The computation of the i so spin-asymmetric nuclear matter is possible with a trivial 
modification of the trial wave function used to project out the ground state of the system. 
The asymmetry is defined by 

where A^ is the number of neutrons and Z of protons. We report in Fig. |4] preliminary 
results of the EOS of nuclear matter for different values of a, using a simplified Hamil- 
tonian containing the AV6' NN interaction. The various curves of the figure from the top 
to the bottom are in the same order of the legend, where it is also indicated the number 
of neutrons and protons for each simulation. 

As displayed in the legend of Fig. IH some calculations were performed using a very 
small number of protons. These results could suffer important finite size effects. The 



study of finite size effect particularly due to the kinetic energy is in progress. It is 
possible to reduce the dependence of the kinetic energy on the number of particles by 
using different kinds of periodic boundary conditions: e.g. the so called twist average 
boundary conditions (TABC)|l25|]. The computation of several EOS by using TABC to 
see the effect to the kinetic energy is in progress. 

CONCLUSIONS AND PERSPECTIVES 

We have briefly shown some recent results obtained using AFDMC theory. This re- 
cently developed Quantum Monte Carlo provides results for the binding energy of light 
nucleirtsl] and neutron drops[l 1] which are is very good agreement with accurate GFMC 
calculations. However, unlike those methods AFDMC is well suited to deal with nucle- 
onic systems with up to hundred of nucleons in the Hamiltonian. 

We have presented AFDMC results for the EOS of neutron matter in the normal 



phase m 311 . In the low-density regime where neutrons form a superfluid phase, we mod- 
ified the trial wave function in order to include important BCS correlations to compute 
the superfluid gap of neutron matter[22]- Preliminary results for the EOS of asymmetric 
nuclear matter have also been presented and discussed. 

It should be stressed that, besides having shown that AFDMC theory opens up the 
possibility of studying the properties of large nucleonic systems, with an accuracy which 
goes much beyond that of other commonly used many-body theories, the results we have 
already obtained indicate serious inadequacies of commonly used nuclear interactions 
in the high density regime, of interest in astronuclear physics. 

The problem of determining the nuclear Hamiltonian, and in particular the effect of 
«-body forces with n > 3 is becoming of primarily importance. We are working on 
developing a new form for the interaction that perturbatively contains the excitation 
of nucleons. The corresponding potential naturally generates many -body forces, that 
should be small in nuclei, but of primary importance in nuclear and neutron matter. 
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